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The gravitational radiation from point particle binaries is computed at the third post-Newtonian 
(3PN) approximation of general relativity. Three previously introduced ambiguity parameters, com- 
ing from the Hadamard self-field regularization of the 3PN source-type mass quadrupole moment, 
are consistently determined by means of dimensional regularization, and proved to have the values 
^ = —9871/9240, K = and C, — —7/33. These results complete the derivation of the general rela- 
tivistic prediction for compact binary inspiral up to 3.5PN order, and should be of use for searching 
and deciphering the signals in the current network of gravitational wave detectors. 

PACS numbers: 04.30.-w, 04.25.-g 



Astrophysical systems known as inspiralling compact 
binaries (ICBs) — two neutron stars or black holes driven 
into coalescence by emission of gravitational radiation — 
are prominent observable sources for the gravitational 
wave observatories LIGO and VIRGO. The appropri- 
ate theoretical description of ICBs is by two structure- 
less point-particles, characterized solely by their masses 
mi or m2 (and possibly their spins), and moving on a 
quasi-circular orbit. Strategies to detect and analyze the 
very weak signals from compact binary inspiral involve 
matched filtering of a set of accurate theoretical template 
waveforms against the output of the detectors. Several 
analyses [ll y, |j| have shown that, in order to get suffi- 
ciently accurate theoretical templates, one must include 
high-order post-Newtonian effects, up to the third post- 
Newtonian approximation (3PN, or ~ 1/c^, where c de- 
notes the speed of light), or even better the 3.5PN ~ 1/c^ 
one. To date, the templates have been completed through 
2.5PN order (for both the phase 4] and amplitude [3]), 
and the specific effects of gravitational wave tails, and 
tails generated by tails themselves, have been added up 
to the 3.5PN order 0. 

Up to now the 3PN-accurate radiation field has only 
been incompletely determined. Previous work at 3PN 
order showed the appearance of ^^ambiguity parameters" , 
due to an incompleteness of the Hadamard regulariza- 
tion (HR) employed for curing the infinite self field of 
point particles. By ambiguity parameter we mean an ar- 
bitrary dimensionless coefficient whose value cannot be 
fixed within HR. In the binary's 3PN Arnowitt-Deser- 
Misner Hamiltonian |3| there initially appeared two am- 
biguities, the "kinetic" ambiguity Uk and the "static" 
one LOg, while in the 3PN equations of motion in har- 
monic coordinates Q there appeared a single ambiguity 
parameter A turning out to be equivalent to w^. The ki- 
netic ambiguity could be resolved by irnposing the global 
Poincare invariance of the formalism [S|, Isf . The ADM 
Hamiltonian and harmonic-coordinates equations of mo- 



tion have been shown to yield completely equivalent re- 
sults [liJ, luj . More recent work using dimensional reg- 
ularization (DR) filially determined the static ambiguity 
to the value ujs = [i3| or, equivalently, A = -1987/3080 
[l3|. The same result was also obtained by means of a 
surface-integral approach [lj|. 

We are concerned here with the problem of the binary's 
3PN radiation field (beyond the Newtonian quadrupole 
formalism), for which three ambiguity parameters, ^, k, 
^, have been shown to appear, coming from the HR of the 
source- type mass quadrupole moment /y of point particle 
binaries at 3PN order [l5j. The terms corresponding to 
these ambiguities are given as follows (see Eq. (10.25) in 
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Here, Gm is Newton's constant, the factor 1/c^ indicates 
the 3PN approximation, and yi, Vi, ai denote the first 
particle's position, velocity and acceleration. The symbol 
1 <-^ 2 means the same terms but with all particles' labels 
exchanged; the brackets () surrounding indices refer to 
the symmetric-trace-free (STF) projection. 

In this Letter we present, for the first time, the values 
of the parameters ^, k and C, and we outline their deriva- 
tion using DR (our detailed investigation will be reported 
in separate papers .16i!l3ill3)- The main strategy is to 
express both the HR and DR results in terms of their 
"core" part, obtained by applying the so-called "pure- 
Hadamard-Schwartz" (pHS) regularization. (Following 
the definition of [13I, the pHS regularization is a spe- 
cific, minimal Hadamard-type regularization of integrals, 
together with a minimal treatment of "contact" ambigu- 
ities, and the use of Schwartz distributional derivatives.) 
The first step of our calculation 1h| is to relate the final 



HR 3PN quadrupole moment, for general orbits, to its 
pHS part: 
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Here, the left-hand side (L.H.S.) denotes the (non- 
circular generalization of the) result of [la|, while the 
right-hand side (R.H.S.) contains both the "core" pHS 
quadrupole moment, and the effect of adding the am- 
biguities (with some numerical shifts coming from the 
difference between the hybrid Hadamard-type regular- 
ization scheme used in []j| and the pHS one). The pHS 
part is free of ambiguity parameters but depends on the 
specific regularization length scales r[ and r'2 introduced 
in the harmonic-coordinates equations of motion |^ . 

The next step is to derive the multipole moments of an 
isolated (slowly-moving) source in d spatial dimensions 
in order to apply DR [l^]. The Einstein field equations 
in d -I- 1 space-time dimensions are "relaxed" by means 
of the condition of harmonic coordinates, d^h^'^ — 0, 
where the gravitational field variable is defined by h^^'^ = 
^/gg^" — "H^" 1 with g^^" being the inverse and —g the de- 
terminant of the usual covariant metric, and with rj^" the 
Minkowski metric (in Minkowskian coordinates). Then, 
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where D^ denotes the flat space-time d'Alembertian op- 
erator, T^'^ the matter stress-energy tensor, A'^" the ef- 
fective gravitational source term (non-linearly depending 
on hP'^ and its space-time derivatives), and t^'^ the total 
stress-energy pseudo tensor of the matter and gravita- 
tional fields. G is related to the usual three-dimensional 
Newton's constant Gn by G = Gat^q"^, where ^0 de- 
notes an arbitrary length scale. We have obtained the 
mass and current multipole moments {II and Jl) of an 
arbitrary post-Newtonian source, generalizing the three- 
dimensional expressions derived in |l9l | to any d dimen- 
sions. The moments II and J^ parametrize the linearized 
approximation in the multipolar-post-Minkowskian met- 
ric exterior to an isolated source 20] . In the case of the 
mass-type moments we find 
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where L = ii- ■ -ig is a multi-index composed of i spa- 
tial indices (^ > 2 is the multipolar order), x^ is the 



STF part of the product of £ spatial vectors [i.e., Xl = 
STF(a;i^ • • • ccj^)], and the time derivatives are denoted by 
a superscript (n). The integrand in Q) is made out of 
the source densities Sc^ = 2[(d - 2)t°° + T'"']/{d - 1), 
Sd c = t"*^ and Yiab = t"-^ , built from the formal post- 
Newtonian expansion, denoted t'^'^, of the pseudo ten- 
sor T^''. For any of these source densities the sub- 
script [£] denotes the infinite post-Newtonian-type ex- 
pansion 5][^](x,t) = Yl=oo4{\^\lcY^^'''^^\^^t), where 
the coefficients are related to the Eulerian F-function by 
ftfe = r(f + £)/[2^''kl r{^+i + k)]. The expression Q 
involves a regularization factor (|x|/ro)^, where B G C 
and ro is a separate "infra-red" (IR) length scale, and a 
particular process of taking the finite part (FP^), which 
constitutes the appropriate d-dimensional generalization 
of the finite part process used in [19| to treat the IR diver- 
gencies of the PN-expanded multipole moments (linked 
to the region |x| -^ +00). 

Since the source densities E, Sa and Safe depend on 
the post-Newtonian expansion of the metric, h , they 
are obviously to be iterated in a post-Newtonian way in 
order to obtain a useful result. At the 3PN order it is 
convenient to parametrize the moments by means of the 
exphcit retarded potentials V, Va, Wat, Ra and X, intro- 
duced when ci = 3 in [3| and generalized to d dimensions 
in ,13.] . Starting from the matter source densities ac^ = 
2[{d-2)T°° + T''']/{d-l), aac = T''^ and aab = T""^ , we 
first define the "linear" potentials V = D^^ [— 47rGcr] and 
Va = D]^^[—ATTGaa], whcrc D]^^ is the usual flat space- 
time retarded operator. The linear potentials are then 
used to construct higher "non-linear" potentials, such as 
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The retardations in the latter potentials are systemati- 
cally expanded to the required PN order. At Newtonian 
order the expression Q reduces to the standard result 
II = Jd'^xpxL + Oic-^) with p = T°°/c^. 

We have used Eq. Q to compute the difference be- 
tween the DR result and the pHS one [l8j . As in the work 
on equations of motion [l^ |l3| , we find that the ambigu- 
ities arise solely from the terms in the integration regions 
near the particles (ri = |x— yi| ^ or r2 = |x— y2| -^ 0) 
that give rise to poles ex l/e (where e = d — 3), corre- 
sponding in 3 dimensions to logarithmic ultraviolet (UV) 
divergencies. We have verified (thanks to the appropriate 
definition of the finite part process FPs) that the IR re- 
gion (|x| -^ -l-oo) did not contribute to the difference DR 
— pHS. The "compact-support" terms in the integrand 
of Q), i.e., the terms proportional to the matter source 
densities a, a a and Gabi were also found not to contribute 
to the difference (thanks to the definition of "contact" 
terms in pHS jl^l)- We are therefore left with evaluating 
the difference linked with the computation of the non- 
compact terms in the local expansion of the integrand in 
Q near the singularities (i.e., ri -^ and r2 -^ 0) that 
produce poles in d dimensions. Let us denote by F^"^^ the 
non-compact part of the integrand of the d-dimensional 



multipole moment Q (including the appropriate multi- 
polar factors such as x^), that is to say we write the non- 
compact part of II as the integral / dflyiF'^'^\ extended, 
say, over two small domains < ri < 7?.i, < r2 < 7^2- 
[At this stage we can set B = in F^'^'^ , and remove the 
FPs prescription.] We write the expansion of F^'^'> when 
ri — > 0, up to any order N ^ in the form (c/. |l3| ') 
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where p, q are relative integers, and ni = (x — yi)/ri. 

In practice the ifp,qS are computed by specializing 
the general expressions of the non-linear retarded poten- 
tials V, Va, Wab, • ■ ■ (vahd for general extended sources) 
to the point particles case in d dimensions. The matter 
source density reads (7 = fii{t) S^'^^x—yi) + l ^^ 2, where 
§{<i) ig Dirac's delta-function in d dimensions, and fli{t) 
denotes a certain function of time, coming from the stan- 
dard prescription for point particles in general relativity, 
and which is computed in an iterative post-Newtonian 
way. The function p,i depends on the potentials V, Va, 
• • • , evaluated at the location of the singular points fol- 
lowing the rules of DR, i.e., by invoking analytic contin- 
uation in d S C. We PN expand the (time-symmetric) 
propagators. For instance, at the IPN order, we use 
D-i = A-i + ■^A-'^d^ + 0{c~^), which yields the solu- 
tion ^ = G/ii I r?"'' -f^^ ^G/2i fcr;'-'*/2c2 (4 - d)] -1-1^ 
2-hC'(c"''), where we denote fc EEr(^)/7r^. Proceed- 
ing further, we insert the previous solution for V into the 
quadratic part of the source term for Wat (whose struc- 
ture is ~ dVdV), expand it when ri — > 0, and then 
integrate term by term in order to find the local expan- 
sion (when ri — > 0) of the corresponding solution, using 
the integration formula A~^ [nf rf] = nf r"'^'^ / [{a - £ + 
2){a + £ + d)] . We generate by this method a particular 
solution of the Poisson-like equation we want to solve, 
and we added the supplementary homogeneous solution 
defined in Ref. 13], whose programs we have re- used for 
the present work. 

The difference VI between the DR evaluation of the 
(d-dimensional) local integral J d'^xF^'^\x), and its cor- 
responding, three-dimensional pHS evaluation, i.e., the 
"partie finie" Pisi,e2 J d'^^F^'^^^H^)j i^ expressible in 
terms of the expansion coefficients of (|5} as (see also [ij]) 
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where the brackets denote the angular average over the 
unit sphere in 3-1- e dimensions (with total volume r22+£) 
centered on the singularity yi. The L.H.S. depends both 
on the regularization length scales si, S2 of the Hadamard 
partie finie, and on the DR regularization characteristics, 
e = d — 3 and £o . 



With this definition, the dimensional regularization of 
the 3PN quadrupole moment (indices L = ij) is ob- 
tained as the sum of its "pure-Hadamard-Schwartz" part, 
and of a "difference" computed according to Eq. ^■. 

^f^^^^^o] = 4f^\si,S2] +I?/„[si,S2;e,4]- At this 
stage a check of the result is that the HR scales si and 
S2 cancel out between the two terms in the R.H.S., so 
that /j ■ depends only on e and £o [its dependence on 
e is of the form of a simple pole followed by a finite part, 
a_i e~^ + oq +0(2)] . Because of this independence from 
•si, S2, we can express our result in terms of the con- 
stants ri and r'2 used as fiducial scales in the final result 
of 0, 0] , see Eq. Q . We can therefore write the DR 
result as 

4f'''\e,e,] = /,J"^)[r;,ri] +I?/,,[r;,r^;£,^o] • (7) 

Let us now impose the physical equivalence between 
the DR result Q and the corresponding HR result ||2Jl. 
In doing this identification, we must remember that 
the "bare particle positions", y^"*™ g^j-^^j y2^''°, entering 
the DR result differ from their Hadamard counterparts 
(used in |8|, ll5|, lla ) by some shifts, which were uniquely 
determined in [l3J, and denoted there ^i{r'i;e,£Q) and 
^2{r2;£,£o) (see Eqs. (1.13) and (6.41)-(6.43) in 13]). 
In the present work, we will denote them by r)i and J72 
in order to avoid any confusion with the ambiguity pa- 
rameter ^. In other words, we impose the equivalence 
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2toi yi''rj(^ -I- 1 <-^ 2 denotes the total 



change induced on the quadrupole moment by the latter 
shifts. We find that the poles ~ 1/e separately present 
in the two terms in the brackets of (|HJ) cancel, so that the 
physical ("dressed") DR quadrupole moment is finite 
and given by the limit shown in |(SJ). 

Finally, by inserting the expressions of the DR and 
HR results given respectively by lO and ^ into Eq. 
(|5|). and by removing the pHS part which is common to 
both results, we obtain a relation for the ambiguity part 
A/ij of the quadrupole moment in ^ in terms of known 
quantities, 

A/„[C+^,K,C+Tfo] = lun[2?/.,[r'i,r^;e,4] 

+ Sr,(r[y^;eJo)Ii3] ■ (9) 

The apparent dependence of the R.H.S. on r'l, rj, s and 
£o, is checked to cancel out. Equation ^ then gives three 
equations for the three unknowns ^, k, C, thereby yielding 
the central result of this work: 



e 



9871 
9240 



, K^Q, C 



T_ 
33' 



(10) 



which finally provides an unambiguous determination of 
the 3PN radiation field by DR. 



We have been able to perform several checks of our 
calculation. First of all, we have also obtained C by con- 
sidering the limiting physical situation where the mass of 
one of the particles is exactly zero (say TO2 — 0), and the 
other particle moves with uniform velocity. Computing 
the quadrupole moment of a boosted Schwarzschild black 
hole at 3PN order and comparing the result with j|. ' in 
the limit m2 = we recover exactly the same value for C 
[iTJ ■ This agreement is a direct verification of the global 
Poincare invariance of the wave generation formalism, 
and a check that DR automatically preserves this invari- 
ance (as it did in the context of the equations of motion 
|l2tll,l |). Secondly, the fact that the value of k is zero has 
been checked by showing that there are no corresponding 
dangerously divergent "diagrams" (in the sense of [21|). 
Finally, we have computed the DR of the mass dipole mo- 
ment, ll , following the same method and using the 
same computer programs, and found that it differs by 
exactly the same shifting of the world-lines, rii{r[;e,io) 
and ?72(?'2;e,^o), as in Eq. ((SJ, i.e., Srjli = mi ri\ + l ^^ 2, 
from the 3PN center-of-mass position associated with the 
equations of motion in harmonic coordinates as given by 
Eq. (4.5) of [13. 

The determination of the values H1U|I completes the 
problem of the general relativistic prediction for the tem- 
plates of ICBs up to 3PN order (and actually up to 3.5PN 
as the corresponding "tail terms" have already been de- 
termined jy]). The relevant combination of the param- 
eters Hl()|l entering the 3PN energy flux in the case of 
circular orbits, namely 9 [l5J, is now fixed to 



^ + 2k + C = -' 



11831 
9240 



(11) 



Numerically, 9 ~ —1.28041. The orbital phase of com- 
pact binaries, in the adiabatic inspiral regime (i.e., evolv- 
ing by radiation reaction) , involves at 3PN order a linear 
combination of 9 and of the equation-of-motion related 
parameter A J22|, which is determined as 



7 1039 

3 ~ 4620 



(12) 



The fact that the numerical value of this parame- 
ter is quite small, 9 ~ 0.22489, indicates, following 
measurement-accuracy analyses [2j , that the 3PN (or bet- 
ter 3.5PN) order should provide an excellent approxima- 
tion for both the on-line search and the subsequent off- 
line analysis of gravitational wave signals from ICBs in 
the LIGO and VIRGO detectors. 
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